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ABSTRACT 

The wide-band Suzaku spectra of the black hole binary GX 339—4, acquired in 2007 February during 
the Very High state, were reanalyzed. Effects of event pileup (significant within ~ 3' of the image 
center) and telemetry saturation of the XIS data were carefully considered. The source was detected 
up to ~ 300 kcV, with an unabsorbed 0.5-200 keV luminosity of 3.8 x 10 38 ergs" 1 at 8 kpc. The 
spectrum can be approximated by a power-law of photon index 2.7, with a mild soft excess and a 
hard X-ray hump. When using the XIS data outside 2' of the image center, the Fe-K line appeared 
extremely broad, suggesting a high black hole spin as already reported by Miller ct al. (2008) based 
on the same Suzaku data and other CCD data. When the XIS data accumulation is further limited to 
> 3' to avoid event pileup, the Fe-K profile becomes narrower, and there appears a marginally better 
solution that suggests the inner disk radius to be 5 — 14 times the gravitational radius (1-sigma), 
though a maximally spinning black hole is still allowed by the data at the 90% confidence level. 
Consistently, the optically-thick accretion disk is inferred to be truncated at a radius 5 — 32 times 
the gravitational radius. Thus, the Suzaku data allow an alternative explanation without invoking a 
rapidly spinning black hole. This inference is further supported by the disk radius measured previously 
in the High/Soft state. 

Subject headings: stars: individual (GX 339-4) — X-rays: binaries 



1. INTRODUCTION 

One of the current issues in black hole (BH) research 
is to measure the spin parameter, a. This is possible if 
we can estimate the radius of the innermost stable orbit 
i?i n , which decreases from 6i? g (i? g being the gravita- 
tional radius) down to 1.235 -R g as a increases from 
to 1. A traditional method of measuring i?j n of black- 
hole binaries (BHBs) is to parameterize the optically- 
thick disk emission component in their X-ray spectra in 
terms of the disk-blackbody (diskBB) model (Mitsuda et 
al. 1984). First applied successfully to GX 339-4 (Mak- 
ishima et al. 1986) and LMC X-3 (Ebisawa et al. 1993), 
this method has been calibrated using BH masses esti- 
mated from companion star kinematics, and confirmed 
to give reliable (e.g., within ~ 30%) estimates of i?i n 
(Makishima et al. 2000) if distance uncertainties can be 
neglected. 

A more modern way to estimate a, applicable also to 
active galactic nuclei, is to utilize iron line profiles, which 
become broadened and skewed due to stronger relativis- 
tic effects as a increases (e.g., Fabian et al. 1989). First 
found from the Seyfert galaxy MCG-6-30-15 by ASCA 
(Tanaka et al. 1995), the broad Fe-K lines were later re- 
ported in a number of Seyfert s using, e.g., XMM-Newton 
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(Nandraetal. 2007). The latest Suzaku results on MCG- 
6-30-15 indicate a high spin parameter of a > 0.917 
(Miniutti et al. 2007). 

Similar high values of a were derived from BHBs 
(Miller 2007, Miller et al. 2009), including GX 339-4 
observed with XMM-Newton and Chandra (Miller ct al. 
2004ab). Analyzing the Suzaku data of this BHB ac- 
quired in 2007 February, Miller et al. (2008), hereafter 
MEA08, confirmed the broad Fe-K feature, and argued 
that the object is an extreme Kerr BH with i?j n ~ R g . 

The wide-band Suzaku spectra of Cyg X-l in the 
Low/Hard state (LHS), in contrast, gave R[ n ~ 15i? g 
(Makishima et al. 2008), via both the disk emission anal- 
ysis and the Fe-K line modeling. Likewise, R[ n ~ 8R 9 
was obtained from the Suzaku data of GRO 1655-40 
(Takahashi et al. 2008). Also, a 2-20 keV Tenma ob- 
servation of GX 339—4 (Makishima et al. 1986) in the 
High/Soft state (HSS), with dominant disk emission, 
gave i?i n = 57 d% km, or 5.6 Q times R g ; here Q = ds/mj, 
d$ is the distance in 8 kpc (Zdziarski et al. 2004), is 
the BH mass in units of a typical value of 7 M@ (Hyncs 
et al. 2004), and R[ n was recalculated applying a correc- 
tion factor of 1.18 (Kubota et al. 1998; Makishima et al. 
2000) and assuming an inclination of i = 30° (Gallo et 
al. 2004). These results imply a ~ 0. 

To examine these apparent discrepancies on a, we 
reanalyzed the same Suzaku data of GX 339-4 as 
MEA03, and found that the X-ray Imaging Spectrome- 
ter (XIS; Koyama et al. 2008) events suffer heavy pileup 
and telemetry saturation. These effects, neglected by 
MEA08, distort the continuum, and indirectly affect the 
Fe-K line shape. 

2. OBSERVATION 

The present Suzaku data of GX 339—4, as used by 
MEA08, were obtained on 2007 February 12, during an 
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outburst which started in late 2006 (Swank ct al. 2006). 
The XIS employed the 1/4 window option and a burst 
option (0.3 sec for XIS0/XIS1, and 0.5 sec for XIS3), to 
achieve a time resolution of 2 sec and a duty cycle of 15% 
or 25%. The Hard X-ray Detector (HXD; Kokubun ct 
al. 2008) was operated in the standard mode. The data 
processing and reduction were performed in the same 
way as MEA08, using the Suzaku pipeline processing 
ver. 2.0.6.13. 

In spite of the 1/4 window and burst options, the XIS 
events piled up significantly (§ 13. 1|) as the source was very 
bright (~ 0.6 Crab in the 2-10 keV Suzaku band). In 
fact, the unscreened XIS event file contains an unusually 
high fraction of Grade 1 events, which are produced when 
pileup occurs. 7 In addition, a variable fraction of the 
CCD frame was often lost due to telemetry saturation; 
this affects the absolute XIS flux. Out of the 10.3 ks of 
XIS0 exposure, only 2.84 ks was free from this problem. 

The HXD-PIN and HXD-GSO spectra, acquired for a 
net exposure of 87.9 ks, were corrected for small dead 
times, but no other correction due to the source bright- 
ness was necessary. We subtracted modeled non X-ray 
backgrounds (NXBs; Fukazawa et al. 2009). The cos- 
mic X-ray background, < 1% of the total counts, was 
ignored. By analyzing the HXD data acquired during 
Earth occultations, we confirmed the NXB models to re- 
produce the PIN and GSO data to within 1%, a typical 
accuracy level (Fukazawa et al. 2009). Since the signal 
becomes comparable to this uncertainty at ~ 300 keV, 
we quote the source detection up to ~ 300 keV. 

In the present Letter, we use data from XIS0, HXD- 
PIN, and HXD-GSO. We do not use data from XIS1 
or XIS3, which are affected more by pileup than XIS0. 
Unless otherwise stated, errors refer to 90% confidence 
limits. 

3. DATA ANALYSIS AND RESULTS 
3.1. Effects of event pileup and telemetry saturation 

To examine in detail the XIS0 data for pileup effects, 
we produced a radial count-rate profile from a 0.5-10.0 
keV XIS0 image (discarding telemetry-saturated frames), 
and compared it with that of a pileup-frce point source, 
i.e., MCG-6-30-15 observed on 2005 August 17. The pro- 
file ratio between GX 339—4 and MCG-6-30-15 decreases 
significantly within a radius of r = 2' of the image cen- 
troid, reaching at the center ~ 25% of the ratio at r > 2'. 
Therefore, the XIS0 data are affected by pileup at least 
within r < 2', and possibly up to 3'. 

Figure [1] shows XIS0 spectra of GX 339—4, accumu- 
lated over different annuli around the source, and divided 
by that outside 3' to visualize shape changes. Pileup 
thus affects the 1-10 keV continuum shape, and produces 
line-like features at 1.8 and 2.2 keV. Although MEA08 
attributed the line features to response uncertainties, it 
fails to explain the fact that the line strength increases 
inwards. These features, appearing at the Si-K edge in 
the XIS and the Au M-edge in the X-ray Telescope (Ser- 
lcmitsos et al. 2007), arc due to rapid changes of the 
instrumental response coupled with pileup. 

Figure [2] shows raw XIS0 spectra accumulated from dif- 
ferent annular regions of the image, including telemetry- 

7 http:/ /www. astro. isas.ac.jp/suzaku/analysis/xis/pilcup/ 
HowToCheckPileup_vl.pdf. 




Energy (keV) 



Fig. 1— The XIS0 spectra of GX 339—4 taken from different 
annuli around the source, 0' — 7.5' (black), 1' — 7.5' (cyan), and 
2' — 7.5' (blue), all divided by that over 3' — 7.5'. The ratios are 
rcnormalizcd to unity at 5.0 keV. The XIS background, though 
inclusive, is negligible. 

saturated frames, and corrected for neither pileup nor 
live-time fraction. They are shown divided by a power- 
law (PL) prediction, which is calculated using respective 
ARFs (Ishisaki ct al. 2007) that properly consider the 
fractional photons falling therein. In all cases, the PL 
was chosen to have T = 2.6 (which approximates the av- 
erage 2-50 keV spectral slope), and a common normal- 
ization of 5 ph cm -2 keV -1 s , absorbed by a column of 
Ah = 5.0 x 10 21 cm~ 2 . Towards the image ccntroid, the 
continua are severely distorted by event pileup. In addi- 
tion, the spectral ratio decreases towards outer regions, 
which suffer increasingly from the telemetry saturation. 
The live-time fractions of XIS0 are 85%, 79%, 60%, and 
43%, from the inner to outer annuli. 

In Figure H the NXB-subtracted HXD-PIN spectrum 
is shown in black. It is further repeated four times, but 
scaled by the XIS live-time fractions so that a pair of 
XIS and HXD-PIN spectra with the same color should 
match. MEA08 used the grey XIS0 spectrum and the 
black HXD-PIN data, but even putting aside the pileup 
distortions, this is an incorrect combination. In MEA08, 
the XIS and HXD-PIN data points match rather acciden- 
tally, because pileup (which increases the highest XIS 
spectral end) and the telemetry saturation (which re- 
duces the XIS normalization) have opposite effects. 

3.2. Analysis of the r = 0' — 4' spectra 

We tentatively fitted the raw XIS0 spectrum from 0' — 
4' (grey one in Figurc[2]), with a PL plus diskBB model, 
but ignoring the 4-7 keV range as MEA08 did. When 
using the PL photon index of T = 2.2 which MEA08 
(supplemented by Miller 2009) found for XIS data, the 
diskBB temperature was T in = 0.8 keV (shown in Figure 
3 (a) in red). The ratio plot (shown in Figure 3 (b) also 
in red) reveals the broad Fe-K line feature reported by 
MEA08. However, the feature becomes much narrower 
and weaker (blue crosses), if we employ T = 2.44 which is 
close to the value of T = 2.4 used by Miller (2009) for the 
HXD-PIN data. In this case, T in = 0.7 keV. Thus, the 
large Fe-K line width claimed by MEA08 depends on the 
employed continuum slopes, which are also indicated in 
Figure^ In fact, the slope should not differ by Ar ~ 0.2 
between the XIS and PIN (suzaku memo-2008-06). 
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Fig. 2. — The XIS spectra from annuli with different radii, di- 
vided by predictions of a common PL model with T = 2.6; 0' — 7'. 5 
(grey), 2'.0-7'.5 (blue), 3'.0-7'.5 (red), and 4'.0-7'.5' (green). 
They are corrected for neither pileup nor telemetry saturation. The 
13—60 keV HXD-PIN spectrum is shown in black. It is reproduced 
after scaling by live-time fractions of the XIS spectra of the cor- 
responding colors: 85%, 79%, 60%, and 43% for grey, blue, red, 
and green, respectively. The two value of T, used in Figure 3, are 
illustrated. 
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Fig. 3.— (a) The XIS0 spectrum of GX 339-4 in vFv form, 
taken from the 0' — 4' region around the source. It is fitted with 
two different sets of PL+diskBB models (see text), ignoring the 
4-7 keV range, (b) The same spectrum, divided by the two models 
given in (a). 



3.3. Analysis of the r > 2' spectra 

Discarding the r < 2' region where the pileup effects 
are severe, we jointly fit the "blue" pair of XIS0 and 
HXD-PIN spectra in Figure H in the 2.4-9.0 keV and 
13.0-60.0 keV ranges, respectively. Though still piled 
up, this XIS0 spectrum (r = 2' — 7'. 5) retains high signal 
statistics. 

The spectra in Figure [2] exhibit a soft (< 8 keV) and 
a hard (20-40 keV) hump, and an Fe-K emission line, as 
found before (Zdziarski et al. 1998; Ueda et al. 1993). 
We hence employed a model consisting of a diskBB, a 
PL, the associated ionized reflection pexriv (Magdziarz 
& Zdziarski 1994), and a laor (Loar 1991) model for the 
Fe-K line with the emissivity index fixed at q = 3.0 and 
the outer disk radius at 400i? g . The metal abundances 
in pexriv were fixed to solar values except that of iron 
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Fig. 4. — Chi-squares of the joint fits to the XIS0 and HXD- 
PIN spectra, shown as a function of the laor inner radius Rp c . 
Blue and red data points use the spectra of the same color as in 
Figure[2] Filled circles show the absorbed diskbb+pexriv+laor fit, 
with the degree of freedom u given on the ordinate. Open squares 
are derived when diskbb is replaced with compbb, with v decreasing 
by 2. 

{Z Fc ), while the disk inclination was allowed to vary over 
i = 25° — 45°. The "constant" parameter was set at 
Cxis = 0.79 (the live-time fraction), for XIS0 with a tol- 
erance of ±0.05, while at 1.07 (Suzaku memo 2007-11) for 
HXD-PIN. We left free all the other model parameters, 
including the absorbing column Nn. 

The fit was acceptable, and behaved as shown by filled 
blue circles in Figure[5]as a function of the laor inner ra- 
dius, R Fe . We thus find i?p c < 3.6i£ g , which is consistent 
with the conclusion of MEA08. The other parameters 

+0.03 



-0.02 

64+ 02 

U -°*-0.03: 

2.53±g$, and Z Fc = 



are: C X is = 0.78 ±0.02, T = 2.68 ±0.02, T in = 0.73 
keV, the pexriv solid angle fraction n/2ir 
its ionization parameter log£ 
l-77±g;i|. Fixing Z Fe = 1.0 worsened the fit by A X 2 ~ 
30. The laor rest-frame energy was E c = 6.71^053 keV, 

(.1-2) 

and its inclination i = 43° _ 4 . However, the Fc-K 
equivalent width (EW), 2301^0 much exceeds the 
value of ~ 110 eV (George & Fabian 1991) predicted by 
the derived T, Q/2tt, and Z Fc . 

The steep slope of T ~ 2.7 indicates that the source 
was in the Very High state (VHS; Miyamoto et al. 1993), 
wherein the disk emission must be Comptonizcd (Kubota 
& Makishima 2004). We hence replaced diskbb with 
a Comptonized blackbody, compbb; the photon source 
it assumes, a single blackbody, may be reconciled with 
the multi-color diskbb formalism by the free Nr. The 
results of this analysis are given in Figure [4] by open 
blue squares. Thus, the local % 2 minimum at R Fe ~ 
10i? g , which existed in the diskbb modeling, became as 
good as the small-i?Fe solution. This large- R Fc solution 
is characterized by a blackbody temperature of 0.54 ± 



0.01 keV, a Compton optical depth of r 
electron temperature of 9.5^'g keV, E c - 



1.09 



+0.06 



-0.12' 

fi oq+0- 22 

°- zy -0.18 



an 
keV 



42° Again, we find that the Fe-line solutions 



and 

are degenerate, depending on the continuum choice. 

3.4. Analysis of the r > 3' spectra 
We finally analyze the "red" spectra in Figure [2] in 



the same way as in § 13.31 under a constraint of Cxis = 
0.60 ± 0.06. Although this XIS0 spectrum (r = 3' - 7'. 5) 
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Fig. 5.— (a) The XISO (r = 3' -7'. 5), HXD-PIN, and HXD-GSO 
spectra (black crosses), simultaneously fitted with the absorbed 
diskbb+pexriv+laor model (red) and presented in a deconvolved 
vFv form. The PL (blue), reflection (green), diskbb (orange), and 
laor (purple) components are also shown, (b) The fit residuals, 
(c) The same as panel (b), but the laor component is set to 0. 

is estimated to be still weakly (< 10%) piled up, we do 
not correct for this, because no established method is 
available yet. The HXD data are the same as in § 13.31 In 
advance, the red XISO spectrum and that of XIS3 from 
r = 3.5' — 7.5', where the count rates are comparable, 
were confirmed to have a constant ratio within ~ 10%. 
We also confirmed with the MCG-6-30-15 data that XIS 
spectra from these outer regions have a constant (within 
~ 10%) ratios to those from the image core. 

As shown in Figure U in red, the results from the 
diskbb and compbb modeling agree better than for r > 
2, with r = 0.16^0 i6 f° r compbb- Below, we exam- 
ine the diskbb4-pexriv+laor fit (filled red circles). In 
contrast to the "blue" spectra (§ I3.3j) . the data now 
prefer the small- Rp c solution as 5.0 < Rp e /R g < 14 
(with minimum at 8.2) at 68% confidence, although 
i?Fe is unconstrained at 90% level due to poor statis- 
tics. The other parameters arc: T- m = 0.72^QQg keV, 

r = 2.67to;o2, ^/2tt = 0.60 ± 0.02, Z Fc = 1.71±g;i|, 
A H = 0.55+g; 2 ] x 10 22 cm~ 2 , and i ~ 33° (the entire 
range of 25°-45° is allowed at the 90 % confidence). The 
line energy E c = 6.681°' 47 keV is consistent with the de- 
rived range of log£ = 3.65;^ 78 . The Fe-line EW, 104±% 
cV, is also now in agreement with the prediction (~ 100 
eV) from V, CI/2tt, and Zy a . Employing i = 30° and the 
correction factor 1.18 (§ Q}, the diskbb normalization 
yields R in = (5T±{\) d 8 km, or R in /R g = {5.6±{%) Q. 
The results remain largely unchanged when fixing Z-p e at 
1.0. 

As shown in Figure |5fa) , the best-fit model has been 
extended to incorporate the 70-300 keV HXD-GSO data, 
resulting in an acceptable (x 2 /v = 108.7/102) simultane- 
ous fit. The model parameters have remained unchanged 
within their statistical errors. Panels (b) and (c) therein 
show the fit residuals. We tried a simple continuum blur- 
ring, and found no important change, but we defer de- 
tailed analysis to later work. 



4. DISCUSSION 

We reanalyzed the Suzaku XISO and HXD data of GX 
339—4 acquired in the 2007 VHS. When the disk emission 
was modeled with diskbb, the "blue" XISO spectrum, 
using r > 2', gave a small Fe-K inner radius as i?Fe < 
3.6i? g (§ I3.3[) . Although this apparently reconfirms the 
high BH spin by MEA08, the Fe line EW is too large 
(§3.3). In addition, a large- Ry& (~ 10i? g ) solution is also 
allowed when considering disk Comptonization that is 
usually the case in the VHS. Thus, the Fe line shapes 
degenerate, depending on the continuum modeling. 

To further eliminate the XIS pilcup effects, we utilized 
the "red" XISO spectrum from r > 3' (§3.4). Then, the 
diskbb and compbb results became more consistent. Im- 
portantly, the large- i?Fe solution (e.g., 5.0 < R Fe /R & < 
14 at la), are preferred, although the small- Rp e solution 
still remains valid at the 90% confidence level. The im- 
plied Fe-K line EW and the center energy are consistent 
with the reflection parameters (§3.4). 

Assuming isotropic emission, the 0.5-200 keV lumi- 
nosity is 3.8 x 10 38 d| erg s -1 , or ~ Q.Ad\m 7 l times the 
Eddington limit (Leci)- Since the VHS have so far been 
observed in other sources over a typical luminosity range 
of (0.2 - l)L Ed (Kubota & Makishima 2004), this indi- 
cates 0.5 < dl/m 7 < 2.5, or 0.7/^m7 < Q < l.Z/y/rwj. 
Considering extreme BH masses of 3 M@ (m-j = 0.43) 
and 15 M Q (m 7 = 2.1), we then obtain 0.5 < Q < 2. 

The diskbb radius R m /R g = (5.61^) Q, found in 
§ 3.4, is in fact a lower limit, because a significant fraction 
of disk photons will be Comptonizcd into the PL by hot 
electron clouds (Kubota & Makishima 2004). Then, the 
true disk area must be a sum of the diskBB area and that 
of the seed photon source. Since the 0.5-200 keV photon 
number in the PL component is ~ 5 times larger than 
that contained in the 0.5-10 keV diskBB emis sion, the 
estimated radius will increase by a factor of y/l + 5, to 
Rin/Rg = (10— 16) Q. The above estimated uncertainty 
in Q then yields R m /R g = (5 — 32). Since we derived this 

range assuming i = 30° (or vcosl = 0.93), uncertainties 
in i can increase the upper bound, but would not affect 
the lower bound by more than ~ 7%. 

The present Fe-line and diskbb analyses consistently 
suggest Ri n /R g > 5, and hence a < 0.4, in contrast to 
the value of a = 0.89±0.04 by MEA08. If our interpreta- 
tion is correct, GX 339-4 is inferred to be spinning only 
mildly (if at all), like Cyg X-l (§[Q Miller et al. 2009). 
In addition, the consistency between the two methods 
reconfirms that Q is relatively close to unity. 

We must admit that the assumption of a single PL 
continuum down to ~ 1 keV might not be warranted. In 
that sense, the value of Ri n /R g = 5.6 Q (§[!]) measured 
with Tenma is more reliable, since it was obtained in the 
HSS where the spectral modeling is much less ambigu- 
ous, and by a non-CCD instrument that is free from any 
pileup. This, together with our estimates on Q, further 
argues against the large spin parameters. 

So far, large values of a were reported for GX 339—4 
based on Chandra data, and from XMM-Newton plus 
RXTE (Miller et al. 2004ab, and Reis et al. 2008). 
However, these measurements have not been examined 
for systematic effects due to the continuum choice em- 
ployed therein. In addition, the former may be limited 
by the continuum bandwidth, and the latter could be 
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still subject to CCD pilcup. 
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